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(Dated: May 2, 2007) 

Long-lived, heavy particles are predicted in a number of models beyond the standard model of 
particle physics. We present the first direct search for such particles' decays, occurring up to 100 
hours after their production and not synchronized with an accelerator bunch crossing. We apply the 
analysis to the gluino (g), predicted in split supersymmetry, which after hadronization can become 
charged and lose enough momentum through ionization to come to rest in dense particle detectors. 
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Approximately 410 pb _1 of pp collisions at ^/s = 1.96 TeV collected with the DO detector during 
Run II of the Fermilab Tevatron collider are analyzed in search of such "stopped gluinos" decaying 
into a gluon and a neutralino (x?), reconstructed as a jet and missing energy. No excess is observed 
above background, and limits are placed on the (gluino cross section) x (probability to stop) x 
[BR(g— ><?Xi)] as a function of the gluino and Xi masses, for gluino lifetimes from 30 (is - 100 hours. 

PACS numbers: 14.80.Ly, 13.85.Rm, 12.60.Jv, 11.30.Pb, 13.85.-t, 14.80.-j 



Split supersymmetry is a relatively new variant of su- 
persymmetry (SUSY), in which the SUSY scalars are 
heavy compared to the SUSY fermions Due to the 
scalars' high masses, gluino decays are suppressed, and 
the gluino can be long-lived. Other new models, such 
as Gauge-mediated SUSY, can also predict a long-lived 
gluino or other heavy, colored, long-lived particles Q. 
The gluinos hadronize into "R-hadrons" [3j], colorless 
bound states of a gluino and other quarks or gluons. As 
studied in Ref . [4j , some 30% of R-hadrons at the Teva- 
tron can become "stopped gluinos" by becoming charged 
through nuclear interactions, losing all of their momen- 
tum through ionization, and coming to rest in surround- 
ing dense material. We present the first direct search 
for the decays of such particles, with deposited hadronic 
energy not in-time with a pp collision. 

A data sample corresponding to an integrated lumi- 
nosity of 410±25 pb _1 fSj, taken with the DO detec- 
tor [6|] from November 2002 to August 2004, has been 
analyzed to search for stopped gluinos. The DO detector 
has a magnetic central tracking system surrounded by 
a uranium/liquid- argon calorimeter, contained within a 
muon spectrometer. The tracking system, located within 
a 2 T solenoidal magnet, is optimized for pseudorapidi- 
ties \r]\ < 2.5, where r\ = — ln[tan(0/2)], and 9 is the 
polar angle with respect to the proton beam direction 
(z). The calorimeter has a central section (CC) covering 
up to |?7| 1.1, and two end calorimeters (EC) extending 
coverage to |7/| ~ 4.2, all housed in separate cryostats Q. 
The calorimeter is divided into an electromagnetic part 
followed by fine and coarse hadronic sections. Calorime- 
ter cells are arranged in pseudo-projective towers of size 
0.1x0.1 in rj x <fi, where </> is the azimuthal angle. The 
muon system consists of a layer of tracking detectors 
and scintillation trigger counters in front of 1.8 T iron 
toroidal magnets (the A layer), followed by two similar 
layers behind the toroids (the B and C layers), which pro- 
vide muon tracking for |?7| < 2. The luminosity is mea- 
sured using scintillator arrays located in front of the EC 
cryostats, covering 2.7 < |?7| < 4.4. The trigger system 
comprises three levels (LI, L2, and L3), each performing 
an increasingly detailed event reconstruction in order to 
select the events of interest. 

We search for stopped gluinos decaying into a gluon 
and a neutralino, \i- The analysis has slightly reduced 
sensitivity for g^qqxl, which may be a large fraction 
of the decays, depending on the SUSY parameters. The 
gluino lifetime is assumed to be long enough such that 



the decay event is closest in time to an accelerator bunch 
crossing later than the one that produced the gluino. For 
the LI trigger to be live again during the decay even if the 
production event was triggered on, this lifetime must be 
at least 30 /is, due to trigger electronics deadtimc. The 
efficiency for recording the gluino decay is modeled as a 
function of the gluino lifetime, up to 100 hours. When the 
decay occurs during a bunch crossing with no other in- 
elastic pp collision, the signal signature is a largely empty 
event with a single large transverse energy (Et) deposit 
in the calorimeter, reconstructed as a jet and large miss- 
ing transverse energy (fir). 

The trigger for each event requires that neither of the 
luminosity scintillator arrays fired. At least two calorime- 
ter towers of size rj x 0=0.2x0.2 with Et>3 GeV are also 
required at LI. Jets are reconstructed with the Run II 
Improved Legacy Cone Algorithm Q with a cone of ra- 
dius 0.5 in rj x <fi space. A reconstructed jet with £r>15 
GeV is required at L3. Offline, we require exactly one 
jot in the event with E>90 GeV, and no other jets with 
Et>8 GeV. The calorimeter requirements in the trigger 
are nearly 100% efficient for events that pass the 90 GeV 
offline threshold. 

To simulate stopped gluino decays, the PYTHIA § 
event generator is used to produce Z+gluon events, with 
the Z boson forced to decay to neutrinos. Initial-state ra- 
diation is turned off, as are multiple parton interactions. 
The spectator particles coming from the rest of the pp 
interaction, such as the underlying event, are removed 
by removing all far-forward particles with \p z / E\ > 0.95. 
The location of the interaction point is placed inside the 
calorimeter, and events are further weighted such that 
the final decay position distribution is that expected for 
stopped gluinos. The radial location of the gluino when 
it decays depends on the way gluinos lose energy via ion- 
ization and stop in the calorimeters. This calculation was 
performed 0] for a distribution of material similar to that 
of the DO calorimeters and a gluino velocity distribution 
as expected from production at the Tevatron. The rj dis- 
tribution is determined by the fact that gluinos would 
tend to be produced near threshold at the Tevatron, and 
that only slow gluinos would stop. The gluinos are thus 
expected to be distributed proportionally to sin 9. More 
than 75% of gluinos that stop have M<1. Because the 
gluinos are at rest and with their spin randomly oriented 
when they decay, the gluon is emitted in a random di- 
rection. Thus a random 3D rotation is applied to the 
simulated particles. 
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The energy of the gluon, which hadronizes and frag- 
ments into a jet, depends on the gluino and neutralino 
masses: E = (M| — M? )/2Mg. We generate four sam- 
ples of stopped gluinos, containing about 1000 events 
each, using a GEANT-based [13] detector simulation and 
reconstructed using the same algorithms as data. They 
correspond to gluino masses of 200, 300, 400, and 500 
GeV, with a neutralino mass of 90 GeV. These samples 
correspond to generated gluon energies of 80, 137, 190, 
and 242 GeV, respectively. Simulated jets are corrected 
for relative differences between the data and simulation 
jet energy scales. The calorimeter electronics sample the 
shaped ionization signal only once per bunch crossing, at 
the assumed peak of the signal for jets originating from 
a pp interaction, but the gluino decay can occur at any 
time with respect to a bunch crossing. So jet energies 
in the simulation are also corrected (downwards) accord- 
ing to a model of this "out-of-time" calorimeter response. 
The average degradation of energy is 30%, although more 
than half of the jets are not significantly degraded. 

The primary source of background is cosmic muons, 
which are able to fake a gluino signal if they initi- 
ate a high-energy shower within the calorimeter. Hard 
bremsstrahlung is responsible for the majority of the 
showers. These showers tend to be very short, since 
they are electromagnetic in nature and thus have small 
lengths compared to hadronic showers. However, some- 
times a wide, hadronic-like, shower can be created cither 
due to deep- inelastic muon scattering, fluctuations of the 
shower, or detector effects. Cosmic muons can usually be 
identified by the presence of a reconstructed high-energy 
muon. A coincidence of muon hits in the B and C layers 
of the muon system, behind the thick iron toroid mag- 
net, is very strong evidence of a muon. The A layer muon 
hits are often also caused by the signal, due to particles 
escaping the calorimeters, so are difficult to use for back- 
ground rejection. Sometimes the muon is not detected, 
due to detector inefficiencies, being out-of-time with the 
bunch crossing, or the limited acceptance. 

Another source of background events is beam-halo 
muons, or "beam- muons." These are muons, synchro- 
nized with the pp bunch crossings and traveling nearly 
parallel to the beam. Often, one or more muon scintilla- 
tor hits can be associated with the muon, and the muon 
is measured to be within Ai<10 ns of a bunch cross- 
ing. Another feature of the beam-muons is that they are 
nearly all in the plane of the accelerator beam. Beam- 
muon showers are also typically very narrow in 0, causing 
this background to be negligible once wide calorimeter 
showers are required. 

Since the trigger requires no signal in the luminosity 
scintillator arrays, nearly all of the pp beam produced 
backgrounds are eliminated. An exception is diffractive 
events with forward rapidity gaps in both the positive 
and negative r\ regions. Typical pp events have a pri- 
mary vertex (PV) reconstructed from tracks which origi- 



TABLE I: The selections applied, and the number of events 
passing in data and for a simulated signal with A/§=400 GeV 
and M^o=90 GeV. 

Selection Data Events Signal Events 



Total 


7199133 


2000 


Exactly one jet (E T >& GeV) 


3691036 


1678 


Jet |?7|<0.9 


2742353 


1505 


Jet E>90 GeV 


202568 


805 


No PV 


198380 


803 


Data quality 


189781 


772 


Jet r] and <f> widths >0.08 


5994 


410 


Jet n90 >10 


1402 


383 


No muons 


109 


357 



nate near to each other along the beamline, where the pp 
interaction occurred. Dijet events in the same data sam- 
ple are studied to understand the $ T spectrum and PV 
reconstruction efficiency for beam-related backgrounds. 
After requiring no PV to be reconstructed and large $ T 
(implicit from the requirement of a single high-energy 
jet), the pp events are negligible. 

Other sources of physics background considered are 
cosmic neutrons and neutrinos, both of which are found 
to be negligible. Cosmic neutrons would have to pen- 
etrate the thick iron toroid. Those neutrons that did 
reach the calorimeter would shower preferentially in the 
outer layers on the top of the calorimeter, which is not 
observed. 

Finally, since the signal process is rare, we also consider 
occasional fake signals caused by detector readout errors 
or excessive noise. We require the jet to be in |?7|<0.9, 
since the forward regions of the calorimeter are observed 
to have more frequent (yet still rare) problems. Also, 
the gluino signal tends to be concentrated in the central 
detector region. Remaining problems are isolated to a 
specific set of runs, detector region, or both, and such 
events are removed. 

The following criteria are used to select events con- 
taining "wide-showers": jet 77-width and ^-width >0.08 
and jet uqq >10, where ngo is the smallest number of 
calorimeter towers in the jet that make up 90% of the jet 
transverse energy. The reverse criteria define a "narrow- 
shower." Criteria are also defined which select events con- 
taining "no-muon" or a "cosmic-muon." An event con- 
tains no-muon if there are no B-C layer muon segments 
in the event, and no A layer segments with A(f>>1.5 radi- 
ans from the jet direction. Cosmic-muon events have at 
least one B-C layer muon segment with |At|>10 ns from 
the bunch crossing time. A candidate stopped gluino de- 
cay event contains both a wide-shower and no-muon. 

To estimate the number of such wide-shower no-muon 
events expected from cosmic muon background, we use 
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Jet energy (GeV) 

FIG. 1: A comparison of the wide-shower no-muon data 
(points) to the expected background from cosmic muons (solid 
histogram) and a simulated signal (dashed histogram). 



the assumption that the probability not to reconstruct 
a cosmic muon in the muon system is independent of 
whether the muon's shower in the calorimeter is narrow 
or wide. A subset of the narrow-shower data sample is 
defined which is nearly devoid of beam-muons by requir- 
ing a shower out of the accelerator plane. This cosmic- 
muon narrow-shower data subset has a similar 77 distri- 
bution to the wide-shower data, and the rj and <p shower 
width distributions are not altered significantly when re- 
quiring a muon. The probability to not reconstruct the 
muon in this narrow-shower data sample is measured to 
be 0.11±0.01, independent of shower energy. This prob- 
ability is applied to the wide-shower cosmic-muon data 
sample to predict the jet energy spectrum of wide-shower 
no-muon background events, as shown in Fig. [T] The 
data agree with the estimated background from cosmic 
muons. There is no significant excess in any jet energy 
range, and the data has the predicted shape in 77 and </>. 

We search for a signal in jet energy ranges with widths 
chosen from the jet energy resolutions of the simulated 
signal samples. The ranges are from M — a/2 to M + 2a, 
where M is the mean jet energy of the sample and a is 
the sample's jet energy RMS. An asymmetric window is 
chosen since the background is steeply falling with in- 
creasing jet energy. 

To first order, the detection efficiency for the decays 
of the stopped gluino signal events can be estimated 
from the simulation, but some effects arc not modeled. 
There is a loss of efficiency at the trigger level from the 
requirement of neither luminosity scintillator array fir- 
ing. If a minimum bias collision happens to occur dur- 
ing the bunch crossing when the gluino decays, a lumi- 




20 40 60 80 100 °0 20 40 60 80 100 

Gluino Lifetime (hours) Time (hours) 

FIG. 2: Left: The trigger efficiency vs. gluino lifetime. Right: 
The instantaneous luminosity profile used to model the trigger 
efficiency. Dashed lines indicate a 50% chance of the store 
occurring. 



TABLE II: The data, background, signal efficiency (for 
stopped gluinos where <7^<?Xi), and expected and observed 
cross section upper limits (at the 95% C.L.) for each jet en- 
ergy range, for a small gluino lifetime, less than 3 hours. 



Energy (GeV) Data Bgnd. 


Eff.(%) 


Exp. (pb) Obs. (pb) 


92.5-104.6 


30 37±3.7 


f.7±0.34 


2.61 


I.8I 


112.4-156.6 


39 40±4.0 4.9±0.98 


0.94 


0.89 


141. 3-213.0 


34 3I±3.f 


6.8±f.36 


0.56 


0.71 


f68.7-270.6 


32 26±2.6 


7.2±f.44 


0.48 


0.75 



nosity scintillator array may fire. The fraction of the 
time this occurs has been measured using cosmic-muon 
events triggered on a jet-only trigger with high thresh- 
old. The efficiency of the luminosity scintillator array 
trigger requirement, averaged over the data set, is 75%. 
The probability to have minimum bias interactions dur- 
ing a given crossing is Poisson distributed, with a mean 
proportional to the instantaneous luminosity, approxi- 
mately 20e30 cm _1 s _1 on average for this data set. A 
detailed model of the trigger efficiency is made as a func- 
tion of the gluino lifetime, for lifetimes up to 100 hours, 
using the typical Tevatron store luminosity profile as in- 
put (see Fig. [2]). Stores typically last ~24 hours with 
a 50% chance of another store following, 6 hours later. 
The current luminosity at the time of the gluino decay, 
and thus the chance to have an overlapping interaction, 
is accounted for. Another source of inefficiency is that 
the trigger is not live all the time, but only during the 
"live super-bunches," which make up 68% of the total 
run time. 

The uncertainties from all sources which affect the sig- 
nal acceptance are added in quadrature, totaling (20- 
25)%. They include the modeling of the out-of-time 
jet response (12%), the data/simulation jet energy scale 
(9%), the rj and radial distributions of stopped gluinos 
[(7-9)%], other geometrical or kinematic acceptances 
(5%), and trigger efficiency [(5-15)%]. 

Given an observed number of candidate events, an ex- 
pected number of background events, and a signal effi- 
ciency in a certain jet energy range, we can exclude at 
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FIG. 3: Top: The expected and observed upper limits on 
the cross section of stopped gluinos, assuming a 100% BR 
of g^gXi and a small gluino lifetime (<3 hours), for three 
choices of the x? mass: 50, 90 and 200 GeV, from left to 
right. Bottom: The upper limits observed on the cross sec- 
tion of stopped gluinos, for various assumptions of the gluino 
lifetime, for a x? mass of 50 GeV. Also shown are the the- 
oretical stopped gluino cross sections (dashed lines, shaded 
area), from Ref. [4|, for the range of assumed conversion cross 
sections. 



the 95% C.L. a calculated rate of signal events giving 
jets of that energy, taking systematic uncertainties into 
account using a Bayesian approach (see Table HI]) . This 
is a fairly model-independent result, limiting the rate of 
any out-of-time mono-jet signal of a given energy. 

From the relation between the gluino and x? masses 
and the observed jet energy, results can be translated 
from the generated set of signal samples to any other set 
of (Mg,M^o) which would give the same jet energy. We 
can therefore place upper limits on the stopped gluino 
cross section vs. the gluino mass, for an assumed x\ mass, 
assuming a 100% branching fraction for g^gXi- These 
can be compared with the predicted cross sections for 
stopped gluinos (which include its production rate and its 
probability to stop) taken from Ref. Q ■ Three curves are 
drawn to represent the large theory uncertainty, resulting 
from the variation of the neutral to charged R-hadron 
conversion cross section used: 0.3, 3, and 30 mb. Fig. [3] 
(top) shows these upper limits for x? masses of 50, 90, 
and 200 GeV, for a small gluino lifetime, less than 3 
hours. If the gluino lifetime is greater than 3 hours, the 



average efficiency of the trigger degrades because signal 
events are not recorded between accelerator stores, and 
the limits become weaker, as shown in Fig. [3] (bottom). 

This is the first search for exotic, out-of-time hadronic 
energy deposits at a high-energy collider. The results 
from 410 pb _1 of Tevatron data are able to exclude a 
cross section of ~1 pb for gluinos stopping in the DO 
calorimeter and later decaying into a gluon and neu- 
tralino. For a x? mass of 50 GeV, we are able to exclude 
Afg<270 GeV, assuming a 100% branching fraction for 
(ji— S'gXii a gluino lifetime less than 3 hours, and a neutral 
to charged R-hadron conversion cross section of 3 mb. 
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